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NUMERICAL  TECHNIQUES  FOR  SMALL  WATERSHED 

FLOOD  ROUTING 

D.  L.  Brakensiek,  A.  L.  Heath,  and  G.  H.  Comers- 
Soil  and  Water  Conservation  Research  Division,  Agricultural 
Research  Service 


INTRODUCTION 

With  the  use  of  high-speed  computers  in  flood  routing  problems,  the  consideration  of  com- 
putational feasibility  is  still  of  concern,  but  at  a  higher  level.  Higher  level  is  used  in  the  sense 
that  a  more  complete  mathematical  approximation  to  the  hydrodynamical  system  is  utilized.  The 
work  reported  in  this  paper  was  carried  out  on  a  small  computer  (IBM  1620,  40  K  memory).2 
For  such  a  computer,  a  balance  between  the  computer  running  time,  problem  input  approxima- 
tion, and  problem  output  accuracy  requirement  is  a  necessity. 

Two  approaches  have  been  taken  in  solving  flood  routing  problems  on  high-speed  digital 
computers.  One  utilizes  flood  routing  methods  such  as  storage  routing,  based  on  continuity  of 
mass  and  a  relation  between  storage  and  outflow  or  both  inflow  and  outflow.  The  other  approach 
is  to  utilize  more  fully  the  basic  hydrodynamic  equations  for  conservation  of  mass  and  mo- 
mentum. The  latter  is  treated  here.  Earliest  American  efforts  using  these  equations  with  a 
high-speed  computer  were  made  by  Issacson,  Stoker,  and  Troesch  (2).-  Harder-  considers  the 
appropriateness  of  several  flood  routing  methods  based  on  these  equations;  however,  his  pro- 
posed method  appears  to  contain  some  difficulties.  These  are  clarified  here.  In  1934,  Thomas  (6) 
discussed  several  degrees  of  approximations  to  these  equations,  but  at  that  early  time  it  was  not 
possible  or  feasible  to  complete  the  computations  for  a  comparison  of  all  methods.  Brakensiek 
and  Comer  (1)  completed  the  quantitative  aspects  of  Thomas's  comparison. 


FORMULATIONS 

As  demonstrated  by  Lighthill  and  Whitman  (3),  many  practical  situations  allow  the  inertial 
terms  in  the  equation  of  motion  to  be  neglected.  An  observable  justification  of  this  is  that  the 
bulk  of  a  flood  wave  travels  much  slower  than  the  full  dynamic  theory  would  predict.  The  numeri- 
cal experiments  by  Issacson,  et.  al.,  (2)  also  indicated  that  the  dynamic  effect  is  soon  damped  out 
by  frictional  forces.  In  this  report,  in  addition  to  neglecting  inertial  terms,  it  is  assumed  that  the 
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momentum   contribution   from  lateral  inflow  can  be  neglected.  The  resultant  equation  for  con- 
tinuity of  mass  and  the  motion  equation  are  written  as 

%x       dt      H  (1) 

and 

dF  "  S°  "  Sf  (2) 

where 

Q  =  flow  rate,  c.f.s. 

A  =  area  of  flow,  ft.2 

H  =  depth  of  flow,  ft. 

q  =  lateral  inflow,  c.f.s./ft. 

S0  =  bed  slope,  ft. /ft. 

Sf  =  friction  slope,  ft. /ft. 

Friction  slop  is  evaluated  by  Manning's  formula, 

Sf  =  (Q/K)2 

where  K  is  total  channel  conveyance 

2/3 
K  =  1.486  R    A,  c.f.s. 

n 
R  =  hydraulic  radius,  ft. 


COMPUTATION    STABILITY 

Equations  1  and  2  were  utilized  without  further  assumptions  for  solution  by  a  computer  pro- 
gram. However,  to  illustrate  the  need  for  a  certain  algorithm  and  to  clarify  the  work  by  Harder, 
additional  assumptions  are  temporarily  made.  If  we  assume  a  rectangular  channel  of  constant 
width  B,  that  lateral  inflow,  q,  is  a  constant  in  time  and  space,  and  conveyance,  K,  is  a  single- 
valued  function  of  flow  depth,  H,  then  equations  1  and  2  can  be  combined  in  the  following  fashion. 
Differentiate  equation  1  with  respect  to  x,  i.e., 


dx2  dxdt 

and  equation  2  with  respect  to  t,  i.e., 

d2H      =  _   2Q     dQ        2Q2     dK 
dxdt  K2    bt  K3    dt 
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See  footnote  4,  page  3. 
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d2H 
Eliminating  ^ — 5    between  these  two  equations,  we  have 
ox  ot 

^Q  =   /  2QB  \    BQ   +  2Q2B    dK 

ax2    \  k2  /  at      k3    at 

Since  K  is  a  single-valued  function  of  H,  we  write 
aK  _  dK  aH 


at 

dH 

at 

and  from 

equation  1 

we 

obtain 

aH 
at  = 

q 

B  " 

1 
B 

aq 
ax 

/Q_    dK\    aQ+    /q_    dK\ 

\BK    dH/   ax        \BK    dH/ 


Substituting  these  we  arrive  at 

aQ  _    K2    a2Q  - 

at         2QB     ax2       \BK    dH/   ax        \BK    dH/  (3) 

Equation  3  represents  a  second-order  parabolic  equation,  under  the  stated  assumptions,  that  is 
equivalent  to  equations  1  and  2.  It  is  the  same  form  as  that  derived  by  Morikawa  (4). 

Following  a  heuristic  approach  from   Richtmyer  (5)  the  stability  constraint  for  an  explicit 
finite  difference  algorithm  for  the  parabolic  equation  3  requires  that 

(At)  <QB  (Ax)2 

K2  (4) 

This  is  the  correct  form  for  Harder's  equation  26.  Contrary  to  the  inference  there,  equation  4 
is  a  very  stringent  constraint,  even  more  so  than  the  one  for  the  hyperbolic  system  which  results 
when  inertial  effects  are  retained.  It  is  also  important  to  note  that  the  stability  constraint  changes 
with  the  solution  being  obtained,  i.e.,  change  the  Q  and  H.  Thus,  a  stable  Ax  and  At  at  one  step  of 
a  problem  may  become  unstable  later  on  in  the  problem.  The  inequality  represented  by  equation 
4  becomes  rather  unrealistic  for  small  watersheds  whereby,  of  necessity,  Ax  is  in  thousands  of 
feet.  One  such  situation  studied  where 

Q  =  4,000  c.f.s. 
B  =  30  ft. 
K  =  1.2  x  105  c.f.s. 
Ax  =  2,000  ft. 

indicated  that 

At    <  30  seconds. 

Lower  values  of  discharge  would  require  even  smaller  At' s.  It  is  apparent  that  the  routing  method 
represented  by  equations  1  and  2  is  not  realistic  with  respect  to  input  for  small  watershed  routing 
or  consistent  with  a  small  computer  when  solved  by  the  usual  explicit  finite  difference  method.  For 
this  reason,  an  implicit  finite  difference  method  was  utilized.  The  resulting  formulation  is  more 
complicated  but  stability  is  unconditionally  guaranteed.  Hence,  for  a  given  watershed  problem  At 
can  be  based  solely  on  considerations  of  desired  accuracy  and  Ax  selected  on  the  basis  of  the  field 
survey  of  the  channel  system. 


To  illustrate  that  computational  stability  is  necessary,  the  simplest  parabolic  equation  is 
used,  i.e.,  the  classic  heat  flow  equation, 


at 

d2Q 
dx2 

The  simpl 

est  ex 

plicit  finite  difference  for  equation  3a 

is  (refer  to  fig. 

1), 

dt  " 

Q4  -  Q3 

At 

a2Q 

/Qs-  Q3 

Q3  -  QA 

1 

_Qs- 

2Q3  +  Qx 

dx2 

\     Ax 

Ax       / 

Ax 

(Ax)2 

(3a) 


Thus,  the  approximation  to  equation  3a  becomes 

Q4  -  Q3  =  Q5  -  2Q3  +  Ql 
At  (Ax)2 


or 


(Q5-2Q3+QiK 


Q4=  _At_ 

(Ax)^  /    "^3  (3b) 


Since  Q4  is  completely  determined  from  Q-  values  of  the  previous  row,  this  formulation  is  called 
explicit.  Further,  it  is  known  that  the  stability  condition  for  the  simple  heat  flow  equation  when 
solved  by  the  explicit  formulation,  Richtmyer  (5)  is, 

(Ax)2       2 

In  table  1  is  shown  what  happens  if  an  error  occurs — say  from  rounding — at  one  point  if, 
for  example,  X  =  1.  The  actual  solution  is  Q  (x,  t)  =  0.  Of  course,  the  initial  and  boundary  con- 
ditions are  also  identically  zero.  An  error  of  1  is  introduced  at  one  grid  point.  Note  that  the  error 
of  1  is  growing  in  magnitude  at  a  rapid  rate. 

In  table  2,  the  same  problem  is  presented  except  that  X  =  1/2.  It  is  apparent  the  error  of 
1  is  being  propagated;  however,  its  magnitude  is  constantly  being  diminished. 

The  above  simplified  example  indicates  that  without  computational  stability,  errors  of  any 
kind  become  so  magnified  that  the  entire  calculation  becomes  gibberish. 


COMPUTATION  FORMULATION 

The  four-point  rectangular  grid  of  figure  1  was  used  to  construct  the  finite  difference  ap- 
proximation to  equation  1,  whereas  only  the  top  two  points  of  the  grid  were  used  for  equation  2. 
Below  are  the  resulting  expressions  for  equations  1  and  2,  respectively: 


Q4  -   Q2  +  Q3  -   Ql   +  A2  -   Ax  +   A4  - 
2Ax                                     2At 

4 
A3  _  1    2    qj 
4    1 

^  -  i  (*  ■*    -  j  *  ) 

These  were  simplified  as  follows 

'%  -   Q2)  +  A2  +  A4  =  a 


H4    -  H2  +  I      ((Q2/K2)2  +  (Q4/K4)2  J  =  0 


where 


Ax  ^      x  '  (6) 


\=    At/Ax 


4 
a=    A(Qi  -  Q3)   +  Ax  +  A3  +(At/2)  2  q. 

1 


0  =  \   (<So>2  +  <So>4) 


The  link  for  simultaneous  solution  of  equations  5  and  6  with  Q  and  A  as  the  dependent  vari- 
ables are  known  relations  (tabular  form)  between  depth  of  flow  and  area  of  flow  and  known  rela- 
tions (tabular  form)  between  area  of  flow  and  total  section  conveyance. 

At  a  given  instant  of  time  for  a  stream  channel  of  N  sections  there  are,  for  the  implicit 
formulation,  (N  -  1)  equation  5's,  (N  -  1)  equation  6's,  and  2  boundary  conditions.  Thus,  there 
are  2N  equations  with  2N  unknowns.  Formally  then,  the  problem  reduces  to  a  solution  of  2N 
simultaneous  equations  at  each  level  of  time.  In  addition  to  the  favorable  stability  feature, 
the  implicit  formulation  insures  continuity  of  flow  being  satisfied  for  the  entire  stream  system. 
It  is  not  unusual  for  explicit  computations  to  result  in  significant  loss  or  gain  of  flow.  However, 
for  these  plus  features,  the  nonlinearity  of  the  simultaneous  equations  requires  an  iterative 
solution  procedure.  Possibly  these  equations  could  be  linearized,  but  the  purpose  of  the  higher 
order  approximations  is  to  retain  the  nonlinear  aspects  of  the  problem. 

Since  flow  is  required  to  be  subcritical,  the  two  boundary  conditions  are  taken  one  at  each 
end  of  the  stream  system.  The  upstream  boundary  condition  is  an  inflow  hydrograph  and  the 
downstream  boundary  condition  is  a  known  rating  function,  i.e.,  a  relation  between  flow  depth  and 
flow. 

All  problems  and  applications  considered  here  require  that  flow  is  occurring  at  the  down- 
stream section  at  the  end  of  the  first  time  increment.  This  applies  to  all  flood  flows  when  there 
is  base  flow,  as  well  as  all  flood  flows  which  occur  concurrently  with  lateral  inflow.  This  does 
not  apply  to  flood  flows  with  no  lateral  inflow,  i.e.,  floods  traveling  down  a  dry  channel.  With 
some  special  adaptation,  the  program  does  work  with  this  situation.  However,  the  results  ac- 
tually are  not  valid,  since  the  basic  equations  are  not  strictly  applicable  because  the  flow  varies 
rapidly,  rather  than  gradually. 


ITERATION  METHOD 

A  step  iteration  procedure  is  used.  Given  a  solution  at  the  itn  section,  the  solution  of  equa- 
tions 4  and  5  is  found  at  the  (i  -  l)tn  section.  Sequential  solutions  continue  in  an  upstream 
direction  to  the  upstream  boundary  condition.  It  is  significant  that  the  solution  sequence  is  in  the 
upstream  direction.  Initial  efforts  in  a  downstream  direction  were  not  successful.  Boundary 
conditions  are  satisfied  as  follows: 

(1)  An  area  of  flow  is  assumed  at  the  downstream  section. 


(2)  Corresponding  to  this  area  the  corresponding  discharge  is  interpolated  from  the  down- 
stream boundary  condition  (rating  function). 

(3)  In  the  sequential  fashion,  the  equations  are  now  solved  for  all  sections,  depth  and  convey- 
ance values  as  needed  are  interpolated  from  the  appropriate  tables  with  area  as  the 
argument. 

(4)  The  calculated  discharge  at  the  upstream  section  is  compared  with  the  known  inflow  (from 
the  boundary  condition). 

(5)  If  a  specified  tolerance  is  exceeded,  a  new  area  is  calculated  for  step  1  and  the  cycle 
repeated.  If  the  tolerance  is  satisfied,  the  particular  set  of  calculated  areas  (depths) 
and  discharges  are  accepted  and  the  time  is  incremented  by  At.  It  is  helpful  to  consider 
that  two  iteration  cycles  are  utilized;  an  inner  cycle  which  solves  equations  5  and  6  on 
a  section- by- section  basis,  and  an  outer  cycle  which  satisfies  the  boundary  conditions. 

The  details  of  the  compuer  program  are  presented  as  the  FORTRAN  II  listing  shown  in 
appendix  I.  This  particular  program  is  for  a  prismatic  channel.  However,  by  reducing  the  number 
of  allowable  sections  and  modifying  several  statements,  a  program  for  nonprismatic  channels  is 
available.  A  number  of  examples  have  been  run  with  this  program;  however,  it  must  be  em- 
phasized that  some  program  difficulties  are  probably  still  undetected.  At  this  stage  the  program 
is  mainly  a  research  tool.  In  the  following  pages,  some  of  the  example  runs  are  presented  to 
point  out  several  aspects  of  flood  routing  by  the  methods  described  above.  In  appendix  II,  format 
information  is  presented.  In  appendix  III,  an  example  set  of  input  is  shown.  In  appendix  IV,  output 
corresponding  to  appendix  III  input  is  presented.  Appendix  V  gives  a  list  of  changes  needed  to 
convert  the  program  to -nonprismatic  channels.  Experience  with  this  program  is  rather  limited 
at  this  time. 


EXAMPLES 

In  all  four  examples,  a  prismatic,  rectangular  channel  30  feet  wide  and  72,000  feet  long  was 
utilized.  All  time  increments  are  1  hour.  Examples  using  critical  flow  at  the  downstream  boundary 
required  some  short  reach  lengths,  down  to  100  feet,  so  as  to  define  the  water  surface  profile  in 
the  drawn-down  region.  The  advantage  of  the  implicit  formulation  was  quite  apparent  in  these 
cases  as  no  evidence  of  instability  was  detected.  It  will  be  noted  that  the  illustrated  section  rating 
loops  are  not  closed  at  the  top.  A  shorter  time  increment  would  have  provided  closing  points. 

1.  The  first  example  illustrates  the  effect  of  the  downstream  boundary  condition.  Two  bound- 
ary conditions  were  studied,  a  critical  flow  and  a  normal  flow  rating  function.  All  other  inputs 
were  the  same.  Figure  2  illustrates  a  spatial  sequence  of  section  rating  functions.  In  figure  2a 
are  the  downstream  boundary  conditions.  In  2b  the  normal  case  shows  a  loop  starting  to  form,  and 
the  critical  case  is  still  single  valued,  but  it  is  shifting  toward  the  normal  case.  In  2c  the  full 
rating  loop  has  developed  and  both  cases  coincide.  Figure  3  shows  the  outflow  hydrographs  for 
the  two  cases.  The  critical  case  shows  slightly  higher  discharges  on  the  rising  side  and  lower 
discharges  on  the  falling  side.  It  also  appears  that  the  rate  of  decline  of  the  recession  is  more 
rapid  for  the  critical  case. 

2.  The  second  example  illustrates  the  effect  of  the  upstream  boundary  condition.  The  in- 
fluence of  the  upstream  boundary  condition  is  illustrated  in  figure  4.  These  are  section  rating 
functions  at  the  same  section,  but  have  resulted  from  different  shaped  inflow  hydrographs. 
Again,  all  other  inputs  were  the  same.  The  inflow  hydrographs  were  triangular  hydrographs 
which  differed  only  on  the  basis  of  time  of  rise  and  time  of  fall.  In  4a  the  upstream  boundary 
condition  was  a  fast  rise  and  slow  fall  hydrograph.  In  4b  the  upstream  boundary  condition  was 
a  slow  rise  and  fast  fall  hydrograph.  Since  the  reference  normal  flow  curve  is  the  same  in 
each  case,  the  effect  of  the  time  distribution  of  inflow  has  been  to  shift  the  rating  function.  The 
influence  on  the  outflow  hydrograph  is  shown  in  figures  5  and  6.  Most  significant  is  the  lag 
between  the  inflow  and  outflow  peak.  The  fast  rise  inflow  produced  a  lag  of  2  hours,  whereas 
the  slow  rise  inflow  produced  a  lag  of  1  hour.  . 
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3.  The  third  example  illustrates  the  effect  of  tributary  flow  into  the  main  channel  at  a  con- 
fluence. A  confluence  was  simulated  by  two  short  reaches  (10  feet).  A  lateral  inflow  was  intro- 
duced at  the  center  section  of  these  two  reaches  so  that  the  lateral  inflow  (c.f.s./ft.)  times  20 
feet  equaled  the  confluence  inflow  hydrograph  at  each  hour  increment.  Thus,  the  confluence  was 
20  feet  wide  and  flow  entered  the  main  channel  uniformly  over  its  entire  width.  As  before,  a  hydro- 
graph  was  also  introduced  as  the  upstream  boundary  condition.  Figure  7b  presents  the  section 
rating  function  at  the  first  section  above  the  confluence  without  confluence  inflow.  Figure  7a  is 
the  section  rating  function  at  the  same  section  with  confluence  inflow.  The  effect  of  backwater  is 
quite  apparent  by  the  shift  from  the  reference  normal  rating  curve.  A  displacement  is  also  noted 
in  the  relative  position  of  the  rising  and  falling  limb  of  the  rating  function.  This  is  because  the 
slope  of  the  water  surface  differs  in  the  case  of  confluence  flow  from  the  case  of  no  flow.  In 
figure  8  are  shown  some  of  the  depth  profiles  above  the  confluence  with  or  without  confluence 
inflow.  In  terms  of  area  inundation  the  confluence  inflow  is  a  major  factor. 

4.  The  fourth  example  presents  a  comparison  of  bed  slopes.  In  figure  9,  the  outflow  hydro- 
graphs  are  shown  for  a  bed  slope  of  0.1  percent  and  1  percent.  Comparison  is  with  the  inflow 
hydrograph.  The  steeper  slope  is  seen  to  be  approaching  pure  translation.  At  the  1-percent 
slope,  the  test  channel  is  nearly  supporting  critical  flow  so  that  this  program  may  not  be  func- 
tioning completely  satisfactorily.  An  indication  is  that  the  outflow  peak  is  higher  than  the  inflow 
peak. 

The  final  example  utilizes  the  programs  as  modified  for  a  nonprismatic  channel.  Data  for 
the  example  are  for  the  storm  of  August  5-6,  1961,  on  watershed  W-l  at  Fennimore,  Wis.  (7). 
Figure  10  presents  the  comparison  of  the  observed  with  the  computed  hydrograph.  More  detail 
on  the  computer  simulation  of  this  watershed  can  be  found  in  Brakensiek.6 


CONCLUSIONS 

1.  For  small  computers  explicit  numerical  formulations  for  flood  routing  may  not  be  feasible. 
Even  for  large  computers  the  usually  small  time  increment  is  not  realistic  with  respect  to  ac- 
curacy in  defining  input.  Output  is  not  generally  required  for  such  small  time  increments. 

2.  The  implicit  algorithm  presented  in  this  paper  appears  to  be  provisionally  limited  as  a 
research  device  but  is  over  coming  the  objections  met  with  the  explicit  method. 

3.  Limited  field  data  checks  have  been  made  with  this  program;  however,  present  results 
appear  to  be  reasonable. 

4.  The  neglect  of  inertia  terms  in  the  equation  of  motion  for  some  flow  situations  may  be 
significant.  In  the  examples  of  this  paper,  a  critical  flow-outflow  requires  considerable  acceler- 
ation which  is  lumped  into  position  head.  However,  even  as  an  approximation,  it  can  be  quite  ac- 
curate in  many  practical  situations. 
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Table  1. — Error  growth  rate  for  an  explicit  formulation  when  X  -  1 


Column/ 
Row 

(1) 

(2) 

(3)                    (4) 
GRID  POINT 

(5) 

(6) 

0 

(Ax) 

2  (Ax) 

3(Ax) 

4(Ax) 

5(Ax) 

Time 
(1)         o 

xo 

20 

0 

0 

0 

20 

(2)     (At) 

0 

0 

3l 

0 

0 

0 

(3)  2  (At) 

0 

1 

-1 

1 

0 

0 

(4)  3  (At) 

0 

-2 

3 

-2 

1 

0 

(5)  4  (At) 

0 

5 

-7 

6 

-3 

0 

(6)  5  (At) 

0 

4-12 

18 

-16 

9 

0 

(7)  6  (At) 

0 

30 

-46 

43 

-25 

0 

(8)  7  (At) 

0 

-76 

119 

-114 

68 

0 

(9)  8  (At) 

0 

195 

-309 

301 

-182 

0 

The  zero  values  in  columns  (1)  and  (6)  are  boundary  conditions. 
The  zero  values  in  row  (1)  are  initial  conditions. 
An  error  of  1  is  introduced  at  column  (3)  row  (2). 
^  Sample  computation:   Q4  =  Qi  -  Q3  +Q5  =  0  -  5  -  7  =  -12 


Table  2. — Error  growth  rate  for  an  explicit  formulation  when  \=  1/2 


Column/ 
Row 

(1) 

(2) 

(3) 
GRID  POINT 

(4) 

(5) 

(6) 

0 

(Ax) 

2(Ax) 

3(Ax) 

4(Ax) 

5(Ax) 

Time 
(1)         0 

Lo 

20 

0 

0 

0 

xo 

(2)    (At) 

0 

0 

3l 

0 

0 

0 

(3)  2  (At) 

0 

1/2 

0 

1/2 

0 

0 

(4)  3  (At) 

0 

0 

2/4 

0 

1/4 

0 

(5)  4  (At) 

0 

2/8 

0 

3/8 

0 

0 

(6)  5  (At) 

0 

0 

5/16 

0 

3/16 

0 

(7)  6  (At) 

0 

^5/32 

0 

8/32 

0 

0 

(8)  7  (At) 

0 

0 

13/64 

0 

8/64 

0 

(9)  8  (At) 

13/128 

0 

21/128 

0 

0 

The  zero  values  in  columns  0  and  5  are  boundary  conditions. 

2  The  zero  values  in  row  0  are  initial  conditions. 
3 
An  error  of  1  is  introduced. 

4  Sample  computation:   Q4  =1/2(Q1  +Q$)  =1/2(0  +5/10)=  5/32. 
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APPENDIX  I  -  FORTRAN  II  PROGRAM 

The  program  listed  on  the  following  sheets  is  written  in  FORTRAN  II.  Testing  is  being  done 
on  an  IBM  1620  (40  K)  machine  with  indirect  addressing. 

Dimension  statements  required  a  floating-point  precision  of  5  decimal  places  and  a  fixed- 
point  precision  of  4  decimal  places  so  as  to  not  exceed  memory  capacity. 

At  this  time  the  STOP  instruction,  statement  No.  2631,  the  258  line,  is  under  study.  Tenta- 
tive results  suggest  a  Branch  to  Statement  No.  263,  or  Statement  No.  240.  If  transfer  is  to 
Statement  No.  240,  line  240,  then  the  following  two  statements  must  be  added  after  Statement 
No.  241,  line  243,  and  before  line  244. 

K3  =  ITS   +1 
IBRS  =  000 

These  are  for  re-initializing  purposes. 

Following  the  mainline  program  the  three  subroutines  are  listed. 
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Parameter  Card  Format 


APPENDIX  II  -  FORMATS 


27  -   29 


30  -   39 


XXX 


XX.XX 


Column 

Field  and 

implied 

decimals 

xx.xxxx 

Symbol 

1  -   10 

ERRC1 

11  -   13 

XXX 

N2  2 

14  -   16 

XXX 

ITS3 

17  -   19 

XXX 

IER 

20  -   23 

XX.XX 

DELA4 

24  -   26 

XXX 

K3 

Explanation 


IBRS 


ERROR5 


decimal  to  calculate  iteration  cutoff,  ERR  1 

entries  per  table 

number  of  columns 

row  following  return  to  initial  inflow 

internal  iteration  area  correction 

last  column  at  which  initial  flow  is  present: 
K3  =  ITS  with  lateral  inflow  or  base  flow 
condition.  K3  =  1  if  a  dry  channel  condition. 

001-flow  is  not  at  last  column. 
000-flow  is  at  last  column. 

Upstream  boundary  condition  tolerance  in 
c.f.s. 


Select  so  that  cutoff  does  not  exceed  the  decimal  part  of  the  field. 

2  Maximum  of  35  entries  per  table. 

3  Maximum  of  70  (columns  or  sections). 
^  Use  area  table  increment. 

5  Select  so  as  consistent  with  required  accuracy,  5  c.f.s.  has  been  used  in  most  examples. 


24 


Input  Data 

Format 

Column 

Field  and 

implied 

decimals 

1     -     8 

9-16 

17  -  24 

25  -   32 

33  -  40 

XXXXXXXX. 

41   -  48 

49  -   56 

57  -  64 

65  -   72 

73  -   75 

XXX1 

76  -   78 

XXX2 

79  -   80 

XX 

1  Start  first  card  001 

2  000       Area 
111        Dept 

222       Conveyance  (total) 

Bottom  section  rated  discharge 
Initial  discharge    i 


Explanation 


Data 


Sequence  number 

Card  type 

Column  or  row  at  which  table  applies 


333 
444 
555 
561 
562 
666 
777 
888 
999 


Initial  area 
SDisch 
SArea 

Initial  inflow 
Slope 


Initial  conditions 


I 


Generally  initial  conditions 

(lateral  inflow) 

(ft./ft.) 
Reach  lengths  Top  to  bottom,  repeat  last  reach  length  (for  each  time  level). 
Lateral  inflow,  time,  and  upstream  discharge 
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APPENDIX  III-  INPUT  LOADING  SEQUENCE  AND  EXAMPLE  INPUT 

1  -  Loading  Sequence: 

(1)  Object  deck 

(2)  Subroutines 

(3)  FORTRAN  II  subroutines 

(4)  Parameter  Card 

(5)  Area  table 

(6)  Depth  table 

(7)  Conveyance  table 

(8)  Terminal  rated  outflow  table 

(9)  Initial  discharges  (t  =o) 

(10)  Initial  areas  (t  =  o) 

(11)  Final  base  flow  (can  be  initial  values) 

(12)  Final  base  flow  areas  (can  be  initial  values) 

(13)  Initial  lateral  inflow  (t  =o) 

(14)  Slope  at  each  section 

(15)  Reach  lengths 

(16)  Lateral  inflow,  time  increment,  and  inflow  at  each  level  of  time  (one  for  each  time  incre- 
ment for  the  duration  of  flood  event) 


26 


ooo<^<-'<-if\)fNjfvi{r>fr>co-tf'-3-ir>ir»>o>o»o«ovO>or~i*-oooo^<>&>o^(>c>CT>osc>o^o<o^<>o^o^o^o^^ 

H(\lfnH(\jmrH(\J(<l  r-H    CM  H    (\l  HMHCgHNH(\|H(\IHf\JHfgHNHNHMHfMH(M 

i-H    CM    CO  r-H    CM  i-l    (NJ 


CO  o 

sD 

o 

m 

o 

o 

in 

o 

in 

o 

i-H 

o 

o 

CM  o 

O 

1—1 

• 

CO 

* 

m 

• 

in 

• 

o 

o 

CO 

43 

f> 

m 

O 

•J- 

r- 

CM 

• 

• 
o> 

i-H 

i-H 

43 

■— i 

vO 

o 

• 

m 

-J 

• 

-J- 

r- 

O 

i— i 

>* 

u 

i-H 

m 

l-H   O 

r- 

c^ 

CM 

X 

o 

f» 

o 

r- 

o 

i-H 

o 

o 

(\ii-i  m 

vC 

tr 

• 

oo 

in 

• 

in 

• 

o 

o 

»-H        (Nl 

*0 

43 

• 

CO 

CM 
CO 

i— < 

• 

•—t 

43 

r-H 

CM 

o 

• 

o 

i-H 

CO 

co 

CM 

C^ 

r- 

CO 

H- 

• 

• 

vO 

o 

r~ 

1— 1 

UJ 

CO 

<J- 

^  o  o 

en 

r- 

CO 

CO 

O 

o 

o 

o 

o 

l-H 

o 

o 

X      o  o 

CO 

43 

CO 

• 

o 

m 

m 

r- 

m 

r» 

o 

o 

O        fM  O 

CO 

• 

• 

r-H 

r» 

tf 

■J" 

• 

i— ( 

• 

i— i 

• 

o 

o 

QC              iH 

CO 

43 

c^ 

lO 

r» 

r- 

CO 

• 

iO 

CM 

CM 

• 

CM 

a 

r— 1 

CO 

nj 

CO 

co 

CM 

<r 

vO 

vO 

a 

• 

CM 

CM 

CO 

CO 

< 

■ 
CO 

CM 

CO 

CO 

QMOO 

r> 

m 

o 

0> 

in 

O 

o 

r- 

o 

!*■ 

o 

r-H 

o 

o 

Q£       in  o 

43 

CO 

CO 

O 

a- 

m 

• 

m. 

• 

o 

o 

<          HJl 

vO 

• 

r-l 

C0 

CM 

0> 

• 

i—i 

CM 

— H 

CM 

o 

i-H 

u 

43 
t-t 

CO 

i— ( 

CO 

o 

O 

co 

• 

•— 1 

iH 

r- 

43 

i£> 

• 

<f 

or 

o 

CM 

0> 

vO 

* 

< 

« 

• 

o 

CO 

a 

CM 

CM 

CM 

f-H  o  o 

CO 

«* 

r- 

f- 

CO 

O 

O 

o 

o 

i-H 

o 

o 

in      (no 

CO 

vO 

i—i 

r-H 

CO 

m 

r- 

in 

r- 

o 

o 

i-l  oo 

en 

• 

• 

C 

vO 

O 

o 

• 

i— i 

• 

r-H 

• 

o 

o 

CO 

O 

m 

CM 

o 

-o 

• 

o 

<n 

CM 

• 

r» 

o 

CM 

i— 1 

CO 

CO 

iO 

in 

«0 

^0 

• 

1— 1 

O 

• 

l-H 

co 
c<^ 

i-H 

00 
CO 
CM 

<r  o  o 

CO 

CO 

CO 

<r 

in 

o 

O 

r» 

o 

t- 

□ 

l-H 

o 

o 

•  o  o 

co 

CO 

CO 

CM 

ifl 

■CO 

o 

m 

• 

m 

• 

o 

o 

■-•  r- 

CO 

en 

• 

i-H 

r- 

CO 

o 

• 

•—I 

CM 

r-l 

CM 

o 

o 

o 

CO 

CO 

CO 

• 

O 

CM 

i-H 

• 

o 

43 

03 

• 

o 

o 

l-H 

CO 

CM 

^H 

in 

CM 

sO 

CO 

i-H 

o 

o 

• 

•-H 

43 

CO 

in 

CM 

• 

PI 

CM  O 

o- 

i—l 

• 

-H 

1— 1 

(MOO 

f- 

CM 

o 

<r 

o 

o 

o 

o 

o 

O 

o 

O    i-H 

l-H 

o 

o 

o 

•  vO  O 

43 

C\J 

m 

CO 

<f 

ON 

m 

m 

r- 

m 

in 

in 

r- 

m 

o 

o 

o 

o 

O 

vO 

CO 

CO 

v£> 

ifl 

vO 

• 

r-H 

i—i 

• 

• 

i-H 

i— i 

• 

• 

o 

o 

o 

r-l 

O 

vO 

• 

<r 

m 

43 

O 

CO 

CM 

r- 

CM 

r» 

• 

• 

CO 

i-H 

o 

CM 

CM 

• 

o 

o 

CM 

43 

CM 

iH 

o 

o 

i-H 

C^ 

o 

<r 

1-1 

• 

o 

• 

0> 

m 

— i 

eg  r-i  o  o 

CO 

r- 

r- 

in 

m 

in 

o 

o 

p- 

43  o 

o 

r^ 

iO 

o 

O    r-i 

r-l 

o 

o 

o 

i-H   •  in  o 

CO 

43 

03 

i-H 

r- 

r- 

i0 

co 

m 

m 

• 

• 

m 

in 

• 

• 

o 

o 

o 

o 

in 

CO 

43 

• 

i—i 

r\j 

i-t 

c- 

CO 

• 

i-H 

r-i 

CM 

IT. 

r-H 

i— i 

CM 

in 

o 

o 

o 

r-H 

r\i 

CO 

43 

43 

i-H 

CO 

CN 

<N 

CM 

o 

-o 

<r 

«o 

<r 

• 

• 

vO 

i-H 

o 
O 

• 

1—1 

• 

0> 

CO 
O 

• 

sO 

a* 

l-H 

i-i 

• 

i-H 

• 

sO 

i— i 

o 

CO 

i-H 

o  o  o  o 

o 

i— ( 

CO 

o 

r- 

o 

o 

O 

f^ 

o 

o 

o 

h- 

o 

o 

O    i-H 

r-H 

o 

o 

o 

O        CO  o 

CO 

CO 

o 

—1 

in 

in 

• 

• 

IT, 

in 

• 

• 

o 

o 

o 

o 

•          <* 

• 

* 

43 

r- 

0> 

i— i 

i-H 

CM 

h- 

r-l 

i-H 

CM 

r- 

o 

o 

o 

CO 

CO 

i-t 

CO 

CM 

• 

43 

m 

i0 

m 

• 

• 

CO 

i-H 

r- 

• 

o 

• 

CO 

r- 

i-H 

m 

CO 

CO 

i-H 

ooooooomo  momo  mooo  m 

Oincor~r»-i^r*-o       r- 

-tvOor-r^r-CM^-cM 

t-H        i-H        cm        co        <»■        <j-         >j- 


ooooooooooooooooo 
ooooooooo 

vOvO\OvO\04343nOvO 

cococococococococo 


ooooooooo 


ooooooooooooooooo 


27 


APPENDIX  IV  -  EXAMPLE  OUTPUT 
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APPENDIX  V-  PROGRAM  CHANGES  FOR  A  NONPRISMATIC  CHANNEL 


(1)  Change  all  dimension  statements: 

(a)  Two-way  arrays  to  (12,  16) 

(b)  One-way  arrays  to  (12) 

(2)  Mainline  Statements 

Line 

028     N  =  1,  delete 

032     K  =  1,  change  to  K  =  ITS 

057^ 

059 

132 

155 

166 

179 

235^ 

(3)  Table  Look-Up  Subroutine 

Line 

012  Change  1  subscript  to  I 


>     Change  N  subscript  to  I 
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